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Abstract 

We consider a 3-parametric linear deformation of the Poisson brackets in classical 
mechanics. This deformation can be thought of as the classical limit of dynamics in so- 
called "quantized spaces". Our main result is a description of the motion of a particle 
in the corresponding Kepler- Coulomb problem. 

1 Introduction 

The set-up of a problem in quantum mechanics usually consists of specifying (a) the structure 
of the algebra of the observables and (b) the Hamiltonian of the problem. If the commuta- 
tion relations in the algebra of the observables are deformations of the usual commutation 
relations between the coordinates and the momenta, the situation is often referred to as 
" dynamics in a quantized space" . Schrodinger was apparently the first to consider such a 
situation: in jj| he calculated the energy levels for the Coulomb potential on a 3-sphere. 
For this problem the term "quantized space" is rather appropriate as the spectrum of the 
distance operator on a 3-sphere is discrete even in the absence of any field. 
One can obtain "classical dynamics in a quantized space" by setting 

Um-[A,B] = {A,B} as h->0, 
h 

where {•, •} are the Poisson brackets of classical dynamics. We prefer to speak about "clas- 
sical dynamics in deformed spaces" (keeping in mind that it is not just the space but the 
whole dynamics being deformed). 

In fact, one does not need a reference to the quantum theory in order to construct a 
deformation of the functional algebra of classical observables. However, we are specifically 
interested in the classical version of the 3-parametric family of deformations considered in 
[pi H U in the context of quantum mechanics. 

This family arises naturally if one considers deformations of the commutation relations 
involving the four space-time coordinates, ten generators of the Poincare algebra (that is, 
infinitesimal translations and proper Lorentz transformations) and the unity element, with 
the property that the resulting algebra contains the Lorentz algebra as a subalgebra. 

In what follows we will consider only three-dimensional non-relativistic problems. The 
usual commutation relations between the dynamic variables are replaced by the following 
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set of brackets satisfying the Jacobi identity: 



{Pi,Xj} = SijI — 



s ' 



{h, Xj} 
{k,Pj} 
{hi lj\ 



£ijk%ki {-^j %i\ — g M 2 ' 
'^ijkPki {-^ \ Pi} ~JjZ "g" 7 
Cijklk, {I, k} = 0. 



{Pi,Pj} = -^, 



* ' 3 J — M 2 ' 




(1) 



Here Xi are the coordinates and pi and U are the components of the momentum and angular 
momentum respectively; the function / is a deformation of the unity function. The three 
parameters L 2 , M 2 and S have the dimensions of area, square of the momentum and action 
respectively. 

In general, these parameters may depend on rotation-invariant functions of the dynamic 
variables such as x 2 ,l 2 , (px) etc. In this case in the resulting non- linear functional algebra 
the Jacobi identities are still satisfied. Here, however, we only consider the case of L 2 ,M 2 
and S being constants. It should be stressed that we do not require the signs of L 2 and M 2 
to be positive. 

In order to obtain the flat space with the usual dynamics on the scale of the Solar system 
one has to assume that the values of L 2 , M 2 and S are huge. In fact, taking M 2 and S to 
be infinite, we get the dynamics in the space of constant curvature equal to 1/L 2 , so L is 
on the scale of the size of the Universe. 

An important feature of the relations ([!]) is their non-invariance with respect to the time 
reversal. Namely, changing the signs of p, I and of the brackets simultaneously one obtains 
relations of the same form but with the sign of the parameter S reversed. 

Various particular cases of the above commutation relations are studied in the literature. 
The case of the space of constant curvature (that is, M 2 , S — 00) has received most attention. 
The first work related to this case is due to Schrodingcr jlj (who assumed the curvature to 
be positive, for the case of negative curvature see ||). Many other works on the subject have 
appeared since; we only mention j|] which contains a discussion of the classical mechanics 
(and, in particular, the Kepler problem) in the space of constant curvature. 

The case S — 00 is due to Yang || . The general case of the above commutation relations 
was introduced in see also [Q. The quantum energy levels for the Kepler-Coulomb 
potential and the harmonic oscillator for this general case are calculated in Q . 

Other approaches to deformed spaces exist: see, for example, ||. 

About the notation: we use bold italic letters for vectors. The scalar and vector products 
of vectors a and b will be denoted by (ab) and [ab] respectively, the absolute value of a is 
denoted by \a\ and its square - by a 2 . 

For simplicity, we always assume the particles to be of unit mass. 

2 The canonical coordinates 
2.1 The canonical coordinates 

Depending on the signs of the constants L 2 , M 2 and L 2 M 2 — S 2 , the functional algebra ([j]) 
is equivalent as a Lie algebra to one of the algebras o(5), o(l,4) or o(2, 3) (||). According 
to the Darboux theorem these algebras can be resolved in terms of four pairs of independent 
canonically conjugate coordinates and momenta and two cyclic variables constructed from 
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the elements of (Q). In the algebraic language, the cyclic variables are the Casimir operators 
of the corresponding algebra. 

The Casimir operators of the algebra under consideration are as follows: 

p 2 2(px) 2 ( 1 1 



K2 1 + L 2 + M 2 S 1 V S 2 L 2 M 2 



and 

K 4 = (Il-[x,p}) 2 



(pi) 2 2( P i)(xi) (xiy 



M 2 S L 2 



The Poisson brackets of K2 and K 4 with any function vanish and their values can be fixed. 
In order to get a correct limit in the usual (that is, nondeformed) case we shall set K2 = 1. 
As for the value of K&, we will only consider the scalar representation of the algebra of the 
observables, that is, set 

II = [asp]. (2) 

This implies that (pi) = (xl) = and, hence, that K4 = 0. The relation (^|) generalises 
the definition of the angular momentum as a function of other variables to the case of the 
deformed space. 

Imposing the relation (||) we reduce the number of Darboux pairs by one. An explicit 
resolution of the algebra (|l]) in terms of three angular variables can be found in the Appendix. 

The equality K 2 = 1 (with I 2 expressed via p, x and I) can be thought of as the equation 
defining the phase space of our dynamical system. In particular, if M 2 , S = 00 and L 2 > 
this is the equation of S 3 x R 3 in R 7 and the relations (|l|) describe the dynamics on the 
3-sphcre of radius L. In fact, the phase space is a manifold whenever S 2 = L 2 M 2 . In 
general, however, it is not smooth: if S 2 7^ L 2 M 2 there is a singularity at / = 0. 



2.2 Problems with spherical symmetry 

For problems with spherical symmetry (such as the Kepler problem considered below) it may 
be of use considering the algebra of observables invariant under rotations. Namely, consider 
the algebra generated by the 5 spherically invariant variables: r 2 = x 2 + p' = p 2 + 
D = (px) + I, I 2 . The commutation relations of the nonlinear algebra generated by 
these 5 variables follow directly from p): 



{/,p' 2 } = 2(^-^), {I,r*} = {I,D} 
{p> 2 ,r 2 } = 4DI, {D,p' 2 }^-2p' 2 I, {D,r 2 } 

The nonlinear functional algebra (H) is generated by 5 elements and has three cyclic 
variables: two of them, namely Ki and come from the initial algebra. The third cyclic 
variable is I 2 which commutes with all the generators of (J3J) . Thus, by the Darboux theorem, 
this algebra can be resolved in terms of one coordinate and one momentum. We choose r 
as the coordinate. (In the non-deformed case this is just the radial variable.) 



T 1 TP 
-- 2r 2 I 



(3) 
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The values of K 2 and K4 being fixed, the quantities D, p' and / depend only on r, 
the corresponding canonically conjugate momentum which we denote by p, and I 2 . Taking 
the squares of the vectors on both sides of (0) and using the fact that K 2 = 1 after some 
manipulations we obtain 

P' 2 = ^(D 2 + l 2 )- (4) 
From the last equality of (||) we get that 




Substituting this expression together with (|J) into the equality K% — 1 we obtain an equation 
on D 

D 2 D 2 + I 2 r 2 2D j2 / 1 1 \ 
1 = — p j ! 1 l 2 I 

r 2 M 2 r 2 L 2 S \S 2 L 2 M 2 ) 
with the general solution 

D = Og. - ^ rHP (| - (1ft. - P) cos JL. (5, 

As a corollary we have 

' - y'(i + ^(i 5 -i4^))( 1 -w) si » b < 6 > 

Thus all elements of the nonlinear functional algebra given by (|^) are expressed in terms of 
the canonically conjugate pair (r, p) and the square of the angular momentum I 2 . 



3 Free motion 

The Hamiltonian of the free motion is chosen so as to commute with infinitesimal translations 
and rotations, that is, with p and I. The simplest such expression which gives the correct 
limit as the parameters tend to infinity is 

It follows immediately from the commutation relations that the vectors p and I are constant. 
Also, 



[xl] [pi] 

L 2 S 

and 



x= Ip -p^- H — — 



• jP_ _ (px) 

S L 2 
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From this 



and, hence, 



■■ 2E 



I = Iq cos 



'2E 



(*-*o) 



Now, taking the vector product of (g) with p one has 

(px)p — p 2 x = I[lp] 

from where 



1? 



LV2E sin 



(t - t Q ) 



P 



\ S P 2 

The value of Jo can be found from the condition K2 = 1 

V 



p 2 



>2E 



(t-t ) \[lp]. 



(7) 



In = 



>2E 



1 + 2E 



L 2 1 
52 ~ M2 



In the case L 2 = 00 the expression for a; should be modified to read 



with 



X = Xq 



[pi] 



S 



n 2 



IoP+ M rr)t+ jrr;pt 



2S' 



1 



AI 2 



P 2(px ) 
S 2 S ' 



The orbits of the free motion are elliptic, parabolic or hyperbolic according to whether 
L 2 is positive, zero, or negative. This should not be surprising given that L~ 2 can be 
interpreted as the curvature of the space. For example, set M 2 , S = 00 and take L 2 > 0. 
Then the elliptic orbits (R) can be obtained from the geodesies on the 3-sphere 



1 = I A 



L 2 



by projecting them orthogonally on a Euclidean 3-space. In particular, the farthest distance 
from the origin reached by a freely moving particle is the same for all orbits and equals to 
L. 
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4 The Kepler problem 



4.1 The equation of the orbit 

The Hamiltonian for the Kepler problem is chosen (as in ||) so as to commute with the 
Riinge-Lenz vector A defined as 

A = [pi] + ™ = j ( P 2 x - (px)p) + ^ (8) 

where r is given by r 2 = x 2 + j-p . 

The commutation relations for the components of A are as follows: 

I 2 

{A l ,A J } = 2(E-—)e ijk l k , 



where E stands for the expression 



1 / , I 2 2al 1 1 \ 



Also, 

{Zj,j4j} = — Ciji-Ak 

and, hence, 

{l 2 ,Aj} = 2e ijk l j A k . 
To see that E commutes with all of the components of A first notice that 

A 2 =a 2 + l 2 (2E~ (10) 

Now, 

2{E,A. l } = {^-,A l }-a 2 {^A l } + ^{l 2 ,A t } = 0. 

Thus E can be chosen to be the Hamiltonian for the Kepler problem. 

The three vectors x, p and A lie in the plane orthogonal to the vector of angular 
momentum. Let us choose the polar coordinates (|a;|,<^) in this plane with the direction 
= being given by the Riinge-Lenz vector. (If A = we take = to be an arbitrary 
ray passing through the origin.) 

Taking the scalar product of (S) with x one has 



A\x\ cos <j) — ar ^ 



I 2 Mh 
Taking the vector product of A and I we have 



(11) 



[Al] = -l 2 p+-[xl]. (12) 



r 



G 



The scalar product of (R2) with x gives 



(px) = -p{[AI\,x)= 1 L ^ . (13) 

Now the orbits can be found as follows. In the expression K2 = 1 one can substitute (px) 
from (|l^) , / — from ( pT| ) and p 2 — from (||) and ( fio| ) . The resulting equation of the orbit 
is as follows: 

2A\x\ sin</> (A\x\ cos 6 - ar) 2 r 2 A 2 - a 2 I 2 , ,, 

H ! — = - — — — I 1 (14) 

IS I 4 L 2 l 2 M 2 S 2 K ' 



In the Cartesian coordinates {x,y) — (\x\ cos0, \x\ sin<p) the equation (14) takes the 
following form: 



2 ,'A 2 + a 2 1 \ 2 fa 2 1 \ 2Ay A 2 2 (\ 



X ( ;4 + T2 I + V ;4 + Tl ) ;c + ;2»^2 ^ 02 



/ 4 L 2 \l 4 L 2 IS l 2 M 2 V S 12 L 2 Af 



2aA / „ „ P , lr , 
= ~ X V y M 2 ' ( ^ 

This equation defines a branch of a plane quartic. The other branch of this quartic 
describes the orbit with the same values of A and I but with the opposite value of a. Notice 
that changing I for —I we obtain, generally, a different orbit. (Recall that the commutation 
relations ([!]) are not invariant under time reversal.) 

Interestingly, it was Giovanni Domenico Cassini (1625-1712) who first suggested that the 
trajectories of planets are quartic curves. Now these curves are known as Cassinian ovals. 
A Cassinian oval is defined as the locus of points with the product of their distances from 
two fixed points equal to a fixed constant. 

4.2 The parametric form for the orbit 

From now on we will assume for simplicity that L 2 is positive. Taking into the account the 
fact that K2 = 1, this implies, in particular, that all orbits are closed. 

The curve ( p"5j ) can be parametrised as follows. The condition L 2 > means that we can 
choose a, b so that 

2 l2 a 2 + A 2 1 

a+b ~ = ^^ + 17 2 



and 



Explicitly, 



Aa 

ab=—. 
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and we take a 2 > b 2 . 

The equation (n5h can be written as 



(ax — br) + vy 



ISv 



Here 



1 I a 2 - A 2 1 



x 2 + A 2 1 V 4A 2 a 



V =2 



I 4 L 2 



I 4 L 2 



I 4 



and 

7 2 = i + r 2 

In the non-deformed limit 7=1 and 



1 1 \ a 



S 2 L 2 M 2 J v 2 



a 2 - A 2 2E 



for E < and v 2 = otherwise. 
We can set 



ax — br — 7 cos 9, 

vy - — - = 7sin6>. 
ISv 



so that 



67 cos + |a| 



\ 



7 2 cos 2 + (a 2 - b 2 ) I I - sin + — — 

ISv 



I 2 

M 2 



1 



a 2 -b 2 



y = -sin 6 + — -j. 



07 cos# + 6— 
a 

.4 



\ 



7 2 cos 2 6»+(a 2 -6 2 ) [I sin 



A 



I 2 \ \ (16) 



ISv 2 J M 2 



In the classical limit the above parametrisation gives the usual trigonometric parametrisation 
of the elliptic trajectory. There exists also a generalisation of the hyperbolic parametrisation. 
Write @ as 



(bx-ar) -[M + Jg^J = s 



with 



1 / a 2 - A 2 1_ 
~T 4 + L 2 



A 2 IV 4A 2 a 2 



I 4 L 2 



I 4 



(17) 
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and 



5 2 = 1 - I 2 
Assume that 5 2 > 0. Then setting 



1 



S 2 L 2 M 2 J n 2 



bx — ar = 6 cosh 9, 
[i>y + = 6smh6. 

we obtain the hyperbolic parametrisation of the orbit. (The case 8 2 < is entirely similar.) 

Notice that ^l 2 v 2 = -Sjp so in the non-deformed limit only one of the parametrisations is 
meaningful for the given values of a and A, the other just giving the equation of the orbit. 
In the deformed case both parametrisations can be meaningful for the same closed orbit. 
In this case the hyperbolic parametrisation is two-valued and the parameter 9 belongs to a 
closed interval. 

In the non-deformed limit the hyperbolic parametrisation remains two-valued with the 
parameter 9 varying over all real numbers. The two branches correspond to the two possible 
values of / = ±1. One branch gives the usual parametrisation of a hyperbolic trajectory. 
The second branch parametrises a trajectory with the same value of A and I but with the 
opposite sign of a. 

4.3 Dependence on time 

Taking the Poisson bracket of x with the Hamiltonian we obtain the equations of motion: 

[ceil \pl] ax ap 



L 2 S Sr M 2 r' 
Substituting / and p from (jl^) and ( |l2| ) respectively, we get 

x=±(Ax-ar)^[xl]-[Al])- [ ^ + ^. 
In particular, for the coordinate y we have: 

V= ^ I I A 2 + a 2 + — ) x - aA 



or, equivalently, 



It follows from (161) that 



V= (ax — br)(bx — ar). 

r 



a\ fb^fcost 



bx — ar a \ WA 
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where A is the expression under the square root in (jig). Thus 



— • M,r Ir)-,/! = ( W(\a\+ lr, ^ i ). (is, 



The function A only depends on sin# so the period of the motion can be found without 
integrating (|lg|): 



T: 2;T '"' 



Z^(a 2 ~b 2 )' 

Notice that for A. I and a fixed, the period only depends on L 2 and not on M 2 and S. 
One can re-write ( |l8| ) as 

/1 ^ 2 + 1 - p (^-r^)(^ 2 -i) 

where 



and /x is defined as in ( |l7| ) so that + v 2 = a 2 — b 2 . Integrating, one gets the following 
parametrisation of the time t by 9: 

M . „ _ „„„. _ _ M 9 + » 8illt - 1 ( jTJ^+fgjg ) . (19) 

" " V\/i + ' 2 (*-ii»)(i-^ 2 )/ 

(The above answer is valid if 

1+ ' J (i-i^) (1 -" 2i2)>0; 

obviuos modifications are needed otherwise in (fL9|).) 

The standard trigonometric parametrisation of t by 9 in the non-deformed Kepler prob- 
lem can be obtained from ([l9|) by setting first S and M 2 to be infinite and then taking the 
limit as L 2 tends to infinity. (Note that as L 2 tends to oo, the parameter \i tends to 0, as 
M=^0 



5 The Hamilton-Jacobi equation 

The spherically symmetric part of the Hamilton-Jacobi equation is obtained by substituting 
(g), H and (|) into the explicit expression for the Hamiltonian function 

1/9 9x 2a a 2 
2E = — {D 2 + l 2 ) + — I- 



r 



10 



Namely, we obtain 



l\/( 1+r2M2 (i-i^))< MV2 - ,2 > si °^ 



where W is the radial part of the action function. One arrives to the same expression 
considering the quasi-classical limit for the Schrodingcr equation in a deformed (quantized) 
space. 

We do not know a direct method of solving this (rather unusual at the first sight) 
equation. In fact, one only needs to know the the second derivative W t e in order to describe 
rasa function of time: r(t) can be found from the equality 



t-tn= / drW, 



'rE- 



Below we obtain an algebraic equation for the second derivative W r E of the radial part of 
the action. 

Let us introduce the following notation: 



1 + r 2 M 2 ( 1 _ 1 ) ) (M 2 r 2 _ P) gin Wr _ a 



^S 2 L 2 M 2 ) ) y ' M M 

In this notation the Hamilton-Jacobi equation becomes 

u> + ^ V = 2Er>-?-^=e. (20) 
Taking the derivative of both sides with respect to E we have 

Finally, taking into the account that sin 2 x + cos 2 i = lwe obtain 

2 2M 2 r 2 - 2 / 2E r 2 - / 1 1 \\ - 

The condition of self-consistensy of the three equations above with respect to only two 
functions u, v leads to an equation for the function F in the form 



(F - OOf = F(F - WO) + {s 2 3 + w 2 3 )F 

17 ' ,! ' S~ 



where s = w = M J and F = F + Asw. 
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6 Conclusion 



The main result of the present paper is a parametric solution of the Kepler problem in a 
deformed space whose construction involves three parameters with dimensions of length, 
momentum and action respectively. 

The orbits of particles in this version of the Kepler problem turn out to be (branches of) 
quartic curves. This should be compared with the analysis of Higgs M who studies the orbits 
in the case where there is only one non-trivial deformation parameter, namely the length L. 
In Higgs' picture the orbits are elliptic, parabolic or hyperbolic just as in the nondeformed 
case. This, however, is due to a special choice of coordinates used in ||. These coordinates, 
though well-suited to the geometry of the problem, have no physical significance. The 
coordinates we use are distinguished among all possible coordinate systems as "the physical 
coordinates". The transformation between Higgs' coordinates and the physical coordinates 
is very simple: in our notations, Higgs' coordinates can be written as Xi/I. 

The fact that the Kepler problem can be solved in the deformed space gives some hope 
that other problems (for example, construction of field theories) could also have reasonable 
solutions in deformed spaces. The discussion of this subject is, however, beyond the scope 
of the present article. 

Appendix: Structure of the algebra (|T|). 

Throughout the section we assume that M 2 > and 



Take the spherical coordinates with the angles (<fi, 9, \) in R 4 - Let (it, U4) be the unit vector 



1 

S 2 



L 2 M 2 



1 



> 0. 



(14,144) = (sin x sin 9 sin 0, sin \ sin 9 cos <fi, sin\ cos 9, cos x) 



and set 



d 



d_ 

86' 



d_ 

~d~4> 



Define the vectors F and L by 



Fi 
F 2 
F, 




and 




T - 9 



respectively. 
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Now, 



/ = (uF) + u 4 p, 

1 . . 

— (p - —as) = pit + [Lu] - u 4 F, 

A ' (23) 



S* 2 L 2 M 
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